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Abstract 



We extract, from a quark model potential that reproduces the number and 
ordering of nonstrange baryonic resonances up to 2.3 GeV, the quantum num- 
bers for the dominant configurations in the ground and first non-radial excited 
states. From the pattern of quantum numbers we identify, from data, spectral 
regularities that allow us to predict the expected high spin low-lying spectrum 
from 2.3 to 3.0 GeV. N — A degeneracies and N parity doublets showing up 
can be interpreted in terms of a simple dynamics. 
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I. INTRODUCTION 



There has been much interest during the last years in the high energy part of the hadronic 
spectrum. The aim is to get a better understanding of the dynamics involved, in particular 
the confinement mechanism in hadrons. On the theoretical side some progress has been 
made. On the one hand, unquenched lattice QCD points out a string breaking in the static 
potential between two quarks [1,2] what should be properly incorporated in the phenomeno- 
logical description of the high energy hadronic spectrum through the coupling to open decay 
channels. On the other hand, the idea of a parity multiplet classification scheme at high 
excitation energies as due to chiral symmetry has been suggested [3] and recently put in 
question [4]. On the experimental side the lack of precise and complete data prevents, at 
the current moment, to establish in a clear-cut way a classification scheme for the highly 
excited hadron spectrum. Hence it may be appropriate to try to extract, from a simple 
dynamical quark model calculation and its comparison to data, some spectral patterns. To 
this purpose we shall consider the light baryon spectrum for which an extensive collection 
of data including high spin state masses, though not very precise in some cases, exists. 

We shall use a nonrelativistic quark potential model. The application of such approach 
for high excited states deserves some comments since the quarks inside the baryons may be 
ultra-relativistic, qq pair creation out of the vacuum becomes more and more relevant as 
the number of open decay channels increases, etc. As it is well-known these shortcomings, 
always implicit in the application of the naive quark model, do not invalidate its usefulness 
to provide information about the pattern of the quantum numbers of the baryon states. 
To get an unambiguous assignment of quantum numbers one should require the model to 
generate, inasmuch as possible, the correct number and ordering of the well established 
experimental states. Different dynamical models providing accurate descriptions of the low 
energy baryon spectrum have been proposed [5] . Most of them rely on non-limited (usually 
linearly dependent on the interquark distance) confinement interactions. As a consequence, 
above 1.9 GeV, there are many more predicted states than observed resonances. This missing 
state problem can be solved by ascribing experimental resonances (in the tiN partial wave 
analysis) to predicted states with a significant coupling to the ttN formation channel [6] . 

Alternatively the missing state problem can be obviated by using a quark-quark screened 
potential as shown in Ref. [7], where a correct prediction of the number and ordering of the 
known N and A resonances, up to 2.3 GcV mass, is obtained. This points out to screening 
as an effective manner to give somewhat account of the coupling to nN channels (lattice 
calculations have unveiled the close relation between screening and the opening of decay 
channels). Such interpretation is supported by the fact that the quantum numbers of the 
dominant configurations (for the ground and first excited states of given J) coming out 
from a screened potential (containing only confinement and a residual one-gluon exchange 
(OGE) interaction) are in perfect agreement with the ones ascribed to data from refined non- 
screened models through the coupling analysis as we have explicitly checked for J < 7/2 
(see Ref. [6] for the identification of the available dominant non-screened configurations) . 

Though one cannot expect, from the minimal screened dynamics employed, to get an 
accurate description of the values of the baryon masses, it is clear its convenience to provide 
an unambiguous assignment of quantum numbers to baryon states beyond the sometimes 
ad hoc arguments that have to be used in non-screened models to identify the experimen- 
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tal resonances when the ttN couphng strengths have similar values for different predicted 

candidates. 

However the smooth screening parametrization employed in Ref . [7] , suggested from old 
lattice data [8], does not keep a direct relation to more recent lattice results [2] pointing 
out a rather abrupt string breaking (rehadronization) transition in the potential. Here we 
deepen this relation. To this purpose the comparison with a sharply saturated potential is 
useful. We analyze the resulting baryon spectrum for both (smooth and sharp) models and 
compare the resulting dominant configurations for any angular momentum. We generate 
patterns of quantum numbers and associate to them spectral regularities observed in the 
data. Prom them we are able to make purely phenomenological detailed predictions for 
resonances above 2.3 GeV, many of them not experimentally identified yet. From data and 
our predictions we identify N — A degeneracies and parity series that have a dynamical 
understanding within our quark model framework. 

These contents are organized as follows. In Sec. H we introduce and discuss the models 
used and their results for the nonstrange baryon spectra. The analysis of the dominant 
contributions and the subsequent derivation of spectral patterns is carried out in Sec. HI. 
Finally, in Sec. IV we summarize our main findings. 



II. THE MODEL 



In the absence of sea quarks lattice QCD predicts the static heavy quark-heavy an- 
tiquark (QQ) interacting potential to rise linearly with the quark-antiquark distance [1]. 
Unquenched (valence + sea quarks) lattice calculations point out the existence of string 
breaking in QCD. This is, the static QQ force becomes screened by intermediate light qq 
pairs, so that when the QQ distance exceeds a critical value (the breaking distance) the QQ 
potential saturates to a constant (saturation energy). In the 1980's a smooth screened po- 
tential was proposed to parametrize this effect [8] . In the last decade such a parametrization 
has been implemented in valence quark models in an attempt to examine the consequences 
of string breaking in heavy quarkonia as well as in the nonstrange baryon spectra [7,9-12]. 

The simplest quark-quark screened potential, containing confinement and one gluon ex- 
change terms, reads (the factor 1/2 has been explicitly written since it reflects the difference 
between the quark-quark and quark-antiquark interactions): 



ar. 



K 

— + 




(1) 



where Vij is the interquark distance, rriij the masses of the constituent quarks, (jjj the spin 
Pauli operators, and Mq is a constant to fix the ground state nucleon mass. The screening 
multiplicative factor appears between parenthesis on the right hand side, n, the screening 
parameter, is the inverse of the saturation distance and its effective value is fitted together 
with the other parameters, a, k, and k^, to the spectrum {fj,~^ ~ 1.4 fm for heavy quarkonia 
and fi~^ ~ 1 fm for nonstrange baryons) [7,12]. 

In the heavy meson sector a good description of bottomonium (the only truly nonrelativis- 
tic quarkonia) is obtained. For nonstrange baryons the model predicts quite approximately 
the number and ordering of the experimental states up to a mass of 2.3 GeV. For the sake 
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of completeness wc show improved results (with more channels entering in the description) 
from Ref. [7] in Figs. 1 and 2. Let us note that the ordering discrepancies, related to the 
relative position of the first excited states of positive and negative parity in N and A, are 
endemic in this kind of quark model treatments; the number discrepancy related to the 
presence of two, instead of one, excited states for N{3/2'^) and N(5/2'^) might have to do 
with the presence of two degenerate excitations at the experimental level of precision. 

More recent lattice QCD calculations [2] show that the QQ potential saturates sharply 
for a breaking distance of the order of 1.25 fm corresponding to a saturation energy of about 
twice the B meson {Qq) mass, indicating that the formation of two heavy-light subsystems 
{B, B) is energetically favored. This information has been implemented in a quark model 
scheme [14] showing that, as a consequence of coupled channels above the physical thresholds 
(corresponding to the opening of decay channels), the description becomes progressively 
less accurate high in the spectrum. Moreover the mixing with the continuum can also 
modify the short-range part of the interaction. Nonetheless an effective (renormalized) 
nonscreened potential continues being useful up to energies not too far above the lowest 
physical threshold. 

These results provide some understanding of the successful application of the smooth 
screened potential of Eq. (1) in Refs. [7] and [12]. For bottomonium, the maximum value 
of the potential, a j \x = 2070 MeV (the h mass has been chosen so that the constant term 
in the potential vanishes), translates into an upper limit for the energy of bound states of 
2m6 + a/iJi = 11400 MeV [12], far above the lowest physical threshold, assuring that the 
potential does not differ much from a nonscreened one (corresponding to Eq. (1) without 
the multiplicative screening factor) for the experimentally known states. As a matter of fact 
one can alternatively use for bottomonium a (linear + coulomb) potential sharply saturated 
at an energy just above the highest resonance experimentally known. 



and maintain the quality of the description [15]. Let us note that this is the closest approach 
one can do to a nonscreened potential if one wants to extend its use high in energy since 
the sharp saturation prevents a proliferation of predicted states without experimental coun- 
terpart. We should then keep in mind that the fitted saturation distance, Tgat-, should not 
be identified with the lattice breaking distance at the physical threshold but as an effective 
parameter to establish the limit of applicability of the model (in fact its fitted value, 1.76 
fm, differs a 30% of the breaking value 1.25 fm). 

In the nonstrange baryon case, the maximum possible value of the energy of bound states 
for the screened potential in Ref. [7], Sm^ + Mq + 3a/ (2/i) ~ 2.4 GeV , just above the highest 
experimental resonance catalogued as four stars by the Particle Data Group (PDG) [13], 
A(ll/2"'") at 2420 MeV, reflects the difficulty of pushing the applicability limit of the model 
beyond this energy. It is interesting to make again a comparison with a sharply saturated 
potential of the form: 




(2) 




fij ^ fsat 



(3) 



where 
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Vsririj) = - avij + 




(4) 



whose parameters, given in Table I, are fitted to get a global description of the nonstrange 
baryon spectrum while choosing Vsat to keep the same value for the saturation energy Sm^ + 
Mq + 'iaTsat/'^ — 2.4 GeV (splitting energy of the three pairs of quarks). We draw this 
sharp potential and the smooth one in Fig. 3. The calculation of the spectrum proceeds 
exactly in the same manner as in Ref. [7], to which we refer for technical details. Let us only 
remind for the moment that two different calculational methods: hyperspherical harmonic 
expansion (HH) and Faddeev, are employed and that convergence of the results is required. 
The resulting spectrum is shown in Figs. 4 and 5. 

It is worth to remark that the presence, in the three-body problem, of two-body thresh- 
olds (for only one quark to be released, values quoted in Table II), apart from the absolute 
three-body one (saturation energy), may represent further constraints in the applicability 
limit of the model to any particular channel. As in Ref. [7] we have also included the 
predicted states close above the thresholds. 

The comparison of Figs. 4 and 5 with Figs. 1 and 2 shows that the sharp potential 
tends quite generally to push upward the highest energy states. However, as we can check 
there is little difference concerning the number and ordering of states. Both models give 
them quite reasonably. In other words the screened potential is quite similar, concerning 
these results, to the closest physical approach to a nonscreened potential, represented by 
the sharp interaction, that takes effectively into account the effect of the baryon decay to 
open channels in order to select the observed resonances. 

Let us notice further that the sharp potential model permits to establish in a direct way 
the subtle connection between the lowest two-body (2bt) and three-body thresholds (3bt): 
if the 3bt is increased (changing for example rgat) and one refits the spin-spin and coulomb 
terms to recover the experimental A—N mass difference, the 2bt is increased, almost linearly, 
in correspondence. Otherwise said, the 2bt (3bt) plus the A — mass difference fixes the 3bt 
(2bt). This is a realization of the idea that the mixing with the continuum (at least partially 
contained in the effective thresholds) affects the whole (in particular the short) range of the 
interaction. Moreover, as the lowest 2bt corresponds to the lowest angular momentum state, 
the experimental presence of a well established high energy-low angular momentum state 
can be used as a lower bound to fix the lowest 2bt. Then this fixing procedure determines 
the maximum energy at which a higher angular momentum state can be predicted within 
the range of applicability of the model (a fine tune of this process has been carried out to 
get the spectrum). 



As shown in Figs. 1, 2, 4 and 5, the predicted nonstrange baryon spectra are quite 
analogous despite the difference between the two potentials. Furthermore the dominant 
configurations entering any ground state and its first nonradial excitation are exactly 
the same for both potentials (this is not the case for higher excitations as discussed below) 
giving for them similar quantitative values that differ at most a 13% from data (within the 
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experimental uncertainties). Hence it makes sense to proceed, from our dynamical models, 
to a quantum number assignemcnt for these states. From the resulting patterns and from 
the experimental values of the masses of the known resonances we shall make predictions for 
higher energy states and we shall justify (as implied by our dynamical model) approximate 
degeneracies and the partial appearance of approximate parity doublets. 



A. ground states 

To express the spatial part of the dominant ground state configurations we shall use 
the hyper spherical harmonic notation, i.e., the quantum numbers {K, L{£, \), Symmetry). 
The so-called great orbital, K, defines the parity of the state, P = {—)^, and its centrifugal 
barrier energy, ^^^(^^ {C — K + ^, p: hyperrradius). £ is the orbital angular momentum of a 
pair of quarks, A stands for the orbital angular momentum of the third quark with respect 

— * — * — * 

to the center of mass of this pair and L {L = i + X) is the total orbital angular momentum. 
Alternatively one can write the parity as P = {—Y~^^, since K — i — X = even > is always 
satisfied. Symmetry specifies the spatial symmetry ([3] : symmetric, [21] : mixed, [111] : 
antisymmetric) which combines to the spin, S, and isospin, T, symmetries {S,T = 3/2 : 
symmetric; S,T — 1/2 : mixed) to have a symmetric wave function (the color part is 
antisymmetric). More precisely, T = 1/2 for and T = 3/2 for A, hence the spatial-spin 
wave function must be mixed for N and symmetric for A. 

It is convenient to group the ground states corresponding to the same dominant config- 
uration in our model as done in Tables III and IV (note that above 2.4 GeV the dominant 
configurations correspond to an extrapolation of the pattern) . 

As a first general rule we have: 

P-H\ (5) 

for the dominant configuration in ground states. This comes, except for A(l/2^), from 
{£ = 0, A) and/or {£, A = 0) components for which L = i + X. For A(l/2+) one has L ~ 
with i = 1 = X instead of £ = = A due to the more attractive spin-spin interaction. 

The second general rule, of general validity for states with J > 5/2, can be formulated 
as follows: 

J > 5/2 : J = L + S with L minimal . (6) 

For states with lower total angular momentum such rule is not always verified. In particular, 
it is not fulfilled by the A^(3/2+) and A^, A(l/2"). For A^(3/2+) the rule would prescribe 
{L = 0,S = 3/2), hence spin symmetric and consequently spatially mixed. Dynamically a 
spatially symmetric {L = 2, S = 1/2) configuration, with a spin-spin attractive contribution, 
is favoured. For A^, A(l/2~), according to the parity expression, L — odd and Eq. (6) 
cannot be satisfied. In all these exceptional cases J = L — S, with S minimal, provides the 
minimum energy. Otherwise said, at the level of dominant configurations, the A^(3/2+) and 
the A^(5/2'^) are degenerate as well as the A^, A(l/2^) and the A^, A(3/2~), respectively. 
These degeneracies are quite well satisfied by data, see Tables 111 and IV. 
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These two rules, complemented with the prescription of having the minimum value of K 
(i.e., the minimum centrifugal barrier) satisfying the symmetry requirements, express the 
minimum energy configuration coming out from our dynamical model. 

It is worth to mention that the same spectral patterns and rules could be also derived 
from a careful look at data in order to apply a general quark model multiplet structure as 
done in Ref. [16]. However, in this case the lack of an underlying dynamics prevents, as 
pointed by the author, a justification for the phenomenological rules derived. 

B. Positive parity states 

For positive parity states the emerging picture consists of two sets (spatially symmetric 
and mixed corresponding to the upper and lower parts of Table III, respectively) of N's and 
A's grouped according to the same dominant configurations. Each A state, except for the 
A(l/2"'"), has a nucleon correspondence. For A(l/2+) its natural nucleon partner should be 
a A^(3/2"'") with (L = 0, S* = 3/2) which, as mentioned above, is dynamically unfavored. 
Actually, our model predicts for iV(3/2"'") (L = 2, 5* = 1/2) a lower mass than for A(l/2+) 
as it is the experimental case. 

It is important to realize that the percentage of the dominant configuration in the sets 
varies: while for the upper part it can be quite different for and A partners (for example 
62% of (2,2, [3]) for A^(5/2+) and 99% for A(7/2+)), for the lower part it is quite approx- 
imately the same and very close to 100% (for instance 99% of (2,2, [21]) for A(5/2+) and 
99% for N{7/2^)). This has to do with the fact that for the lower part series only a spatially 
mixed wave function is allowed whereas for the upper part N series there are, except for 
A^(l/2+), two possibilities open for the same {K, L) values: spatially symmetric (spin-isospin 
symmetric) and spatially mixed (spin-isospin mixed). 

For the upper part set the lowest lying states, iV(l/2+) and A(3/2+), with a dominant 
configuration (0, 0, [3]), have an energy difference of 292 MeV, entirely due in our model to 
the spin-spin potential, attractive for S — 1/2 and repulsive for S — 3/2, since for both 
states the percentage of the dominant configuration is almost the same (98% for A^(l/2'^) 
and 100% for A(3/2+). Let us keep in mind that c/izra/ contributions from meson exchanges 
are taken into account through the effective quark-gluon coupling constant fitted to the 
A(3/2+) — iV(l/2+) mass difference. For the next dominant configuration, containing the 
A^(5/2+) and the A(7/2"'"), whose masses are quite well reproduced by our model, one 
has to consider the configuration mixing effect in A^(5/2+) (62% of (2,2, [3]) and 34% of 
(2,2, [21])). Let us realize though that both configurations would give the same energy for 
an hypercentral (only p— dependent) potential. As our (confinement -|- coulomb) potential is 
quite close to an hypercentral one, the difference between A^ and A partners is again due in 
our model to the spin-spin interaction. We predict 199 MeV. This means a slow variation of 
the spin-spin contribution when increasing K and L for spatially symmetric states (actually 
the 62% of the A(3/2+) - A^(l/2+) mass difference is 181 MeV). The experimental A - A^ 
energy difference is bigger, A(7/2"'") — N{5/2^) ~ 270 MeV, but also seem to vary slowly: 
A(ll/2+) - A^(9/2+) ~ 200 MeV, A(15/2+) - A^(13/2+) ~ 250 MeV. This may be reflecting 
a larger probability for the spatially symmetric component than given by our model (let us 
point in addition that the more cumbersome numerical procedure for J > 9/2 states due to 
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the presence of thresholds, makes less accurate our predictions, indeed the calculated energies 
for A^(9/2+) and A(ll/2+) are inverted respect to data). These quantitative differences 
should be expected given the shortcomings of our dynamical model. However the spectral 
pattern derived is very useful to make manifest a regularity in the data: the experimental 
and A energy steps when going from one dominant configuration to the next one remain 
also constant about 500 MeV: N{9/2+) - N{5/2+) ~ 540 MeV, A(ll/2+) - A(7/2+) ~ 470 
MeV, iV(13/2+) - iV(9/2+) ~ 480 MeV, A(15/2+) - A(ll/2+) ~ 530 MeV (this energy 
step can be associated, from our model, to AK = AL = 2). The extension of this pattern 
to further steps allows us to predict that a A'"(17/2+), if existing identifiable experimental 
resonances at such high energies, would have a mass of 3200 MeV. 

Regarding the lower part set we should notice that the same configuration percentage 
for N's and A's would imply a degeneracy, term by term, if only spatial (neither spin nor 
isospin dependence) interactions were considered, N{7/2~^) : A(5/2+); iV(ll/2+) : A(9/2"'"), 
... As a matter of fact, as for this set the spin-spin interaction (the only non-spatial in 
our model) gives little contribution, due to its short-range character and the spatial mixed 
symmetry involved (with K = L > 2), the difference between the N{7/2~^) and A(5/2'^) 
masses (the only states clearly below threshold) is in our model less than 5% what seems 
to be corroborated by data (this could be an indication that isospin dependent interactions 
do not contribute significantly as well). Concerning the energy step we expect, from our 
analysis above {AK = AL = 2) and the scarce data (A(9/2+) - A(5/2+) ~ 400 MeV) an 
approximate constant value around 400—500 MeV. From it there would appear a iV(ll/2"'") 
at about 2450 MeV (the energy of its two stars A(9/2+) partner is experimentally in the 
interval 2200-2500 MeV), a A(13/2+) and N{15/2+) at about 2900 MeV and a A(17/2+) 
around 3300 MeV. 

All predictions of our model up to 3 GeV, from this subsection and next subsections 
III C and III E, are summed up in Table VII. 

C. Negative parity states 

For negative parity states the results are organized in Table IV. The upper part 
corresponds to spatially mixed configurations for which we expect, from our model, a 
rapid decrement in the spin-spin energy difference between N and A partners when go- 
ing from K — L — ItoK — L>3 and, as a consequence, all the states with the 
same dominant configuration should become close in mass. The data seem to follow 
this prescription: mass differences of 180 MeV at most for K = L = 1, 60 MeV for 
K = L = 3.... Again the spectral pattern makes manifest the energy step (400—500 MeV) 
regularity in the data: (iV(9/2-) - N{5/2-) = 575 MeV; A(7/2-) - A(3/2-) = 500 MeV, 
A^(ll/2-) -iV(7/2-) = 410 MeV) as it corresponds to AK = AL = 2. Therefore one would 
expect a iV(13/2-) and a A(ll/2-) at about 2650 MeV and a A^(15/2-), a A(15/2-) and 
a iV(17/2-) at about 3050 MeV. 

The lower part of Table IV contains an independent series of A's without N's partners. 
All states are above their lowest threshold. Prom the analysis above we can expect again a 
step energy, corresponding to AK — AL — 2, of around 400—500 MeV that turns out to be 
compatible with data within the experimental errors. Prom it a A(17/2^) would be located 
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at approximately 3200 MeV. 



D. N-A approximate degeneracies 

According to our model, the masses of the N{J~) and A(J~) in the left hand side of 
Table IV should become degenerate as the effective spin-spin interaction becomes weaker 
and the potential approaches an hypercentral one. This seems also to be the experimental 
situation. 

On the other hand, the N{J~^) and A{J~^) on the right hand side of Table 111 should 
become degenerate as the spin-spin contribution for 5* = 3/2 (significantly smaller than for 
S = 1/2 in spatially symmetric states) is not quantitatively much different for spatially 
symmetric and mixed states. This is confirmed by the scarce data available (J^ = 7/2+). 

These results, from data and our predictions, are picked up in Table V. We can express 
its content in a schematic simplified manner through: 

An + 3 

Ar(j+'-) ~ A(J+'-) for J= — — with n= 1,2,3... (7) 

Let us realize that for positive parity Eq. (7) corresponds to a rule derived in Ref. [16] from 
a phenomenological mass formula obtained from the assumption of linear Regge trajectories. 



E. Excited states 

For excited states in general the situation is much more cumbersome. Though both mod- 
els, smooth and sharp, provide a similar pattern a more careful analysis of the configurations 
entering the states makes clear significant differences. For instance the fourth A^(l/2+) state 
has for the smooth potential a (2,2, [21]), 5" = 3/2 component degenerate (at the level of 
precision of the model) with the (2,1, [111]), S — 1/2, while the same state corresponds 
uniquely to the last configuration for the sharp potential. For A(3/2"'") in the sharp case, 
the radial excitation of (0, 0, [3]), S = 3/2 and the (2, 2, [21]), S = 1/2 configurations are de- 
generate, whereas in the smooth case, the radial excitation of (0, 0, [3]), 5 = 3/2 is the first 
excited state and it is not degenerate, etc. Moreover, the first radial excitations of A(3/2+) 
at 1600 MeV and iV(l/2"'") at 1440 MeV are badly predicted. In contrast to this inadequacy 
to deal with the radial excitations and the strong model dependence of the configurations 
for higher (second and up) excitations, our dynamics provides a very simple picture for the 
first non-radial excitations that follows the trend of data. The rule, coming from the absence 
of spin-orbit and tensor forces in our model, and again of general validity for J > 5/2, is: 

The first non-radial excitation of N, A( J) and the ground state of N, A( J-l-1) respectively, 
are almost degenerate. (8) 

This rule is satisfied by existing data at the level of 3%, Table VI. From now on we will denote 
by a black dot the first non-radial excitation of any state. The only significant deviation from 
the previous rule, A(5/2-)*(2350) : A(7/2-)(2200) (7% of difference), corresponds to two 
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resonances catalogued as one star by the PDG. Indeed the degeneracy can be accomphshed 
within the experimental uncertainties, see Fig. 5. 

Assuming this approximate degeneracy we can guess, from data and from our phe- 
nomenological ground state predictions above, the masses for the first non-radial excited 
states with a question mark in Table VI. These predictions, up to ~ 3 GeV, are summed 
up in Table VII. 

For J < 5/2 the rule is only satisfied, according to our model, by A^(3/2^) and A(3/2+). 
For some of the remaining states our model gives also definite relations between their first 
non-radial excitations and other ground or excited states. These relations are quite rea- 
sonably satisfied by the data too: A^(l/2-)'(1650) : A^(3/2-)*(1700) : A^(5/2-)(1675), 
A(l/2-)'(1900) : A(3/2-)'(1940), A(l/2+)*(1910) : A(5/2+)*(2000) : A(7/2+)(1950). 

It is worth to remark that some of our findings in Table VII might be in fact already 
accomplished by current data although being masked by the experimental uncertainties (see 
Figs. 4 and 5). For example the A^(7/2+'~)* could be, within theoretical plus experimental 
errors, in the N{7/2+'-), the A(7/2~)* in the A(7/2-) and so on. We should also realize 
that some of the predicted resonances could have small couplings to formation channels what 
would make difficult its detection (see for example Ref. [6] regarding A^(ll/2+), A^(13/2~) 
and N{15/2+)). 

F. Parity series 

A look at Tables III, IV and VII shows some approximate N{J^) parity doublets. Thus, 
for J > 5/2 the N'^ ground state series in the upper part of Table III and the A'"" ground state 
series in the upper part of Table IV appear almost degenerate, term by term, for example 
A^(5/2+)(1680) : A^(5/2-)(1675), as shown on the left hand side of Table VIII. From our 
model we can propose a qualitative systematic for this doubling: the bigger repulsion in 
the positive parity state with respect the negative parity one, due to AK — 1, AL = 1, is 
compensated by the bigger attraction in a spatially symmetric S = 1/2 configuration (i.e., 
AS* = —1). We can quantify these effects from our data analysis above. We can estimate 
an increment of the repulsion energy of ~ 250 MeV (from ~ 500 MeV for AK = 2, AL = 2, 
as previously discussed) and an increment in the attraction, when going from S = 3/2 to 
S — 1/2 of approximately the same amount. So we can predict further parity doublings, 
A^(13/2+) : A^(13/2~), ... Moreover the apphcation of our preceding rules for the first non- 
radial excitations gives rise to an equivalent approximate parity doublet series, right hand 
side of Table VIII. 

For A's (J > 5/2) although an inspection to data in Tables III and IV might also suggest 
some approximate degeneracies: the A( J+'^) ground state series in the lower parts of Table 
III and Table IV could be almost degenerate, term by term: A(5/2+)(1905) : A(5/2-)(1930), 
A(9/2"'")(2300) : A(9/2~)(2400)...., we expect our model not to support this suggestion. 
Indeed going from the positive to the corresponding negative parity state we expect the 
centrifugal repulsion getting increased {AK = 1) and the spin-spin interaction changing 
from attractive {S = 1/2) to repulsive {S = 3/2), hence the J' states should be higher in 
mass than the J"*" ones. 

We should realize that the approximate A^ parity doublets appear discretely, every 500 
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McV. This does not preclude the existence for other energies of additional (A^ and/or A) 
parity doublets involving higher excitations but we cannot carry out for them a trustable 
analysis. 

IV. SUMMARY 

We have performed a phenomenological analysis of the nonstrange baryon spectrum be- 
tween 2.3 and 3 GeV. Our starting point is a quark model potential, containing confinement 
and OGE interactions saturating at 2.4 GcV, that gives the correct number and ordering 
of known resonances up to 2.3 GeV. The study of the dominant configurations entering the 
ground state and the first non-radial excitations provides us with quantum number spectral 
patterns. The extrapolation of these patterns to higher energies allows the prediction, from 
existing data, of ground (J > 11/2) and first non-radial excited (J > 7/2) state masses 
in the quite experimentally uncertain energy region from 2.3 to 3 GcV. It is precisely this 
uncertainty, rcfiecting the transition to the continuum, what makes our picture plausible 
despite the blown up number of new states. 

The quahtative adequacy of our model description to data suggests that the effec- 
tive (confinement -l-coulomb-l-spin-spin) dynamics employed (we have estimated numerically 
small contributions of conventional meson exchange potentials to the mass of J > 5/2 
baryons) contains the essential ingredients to give account of the observed regularities in the 
spectrum. To this respect we should emphasize i) the usefulness of screening as an effective 
mechanism in quark potential models to provide an unambiguous assignment of quantum 
numbers to states (this is related to the connection of screening to decay channel effects as 
shown by lattice calculations) and ii) the importance of using the hypersphcrical harmonic 
basis as the more convenient to deal, from the point of view of the physical interpretation, 
with the quasi-hypercentral potential considered. 

Prom data and our predictions we find on the one hand approximate N — A degeneracies 
and on the other hand series of N parity doublets involving ground or first non-radial 
excitations exclusively. We get from our minimal dynamics some understanding about how 
these degeneracies come out. 

We should finally remark that our results do not rely on any symmetry assumption in 
QCD, but on a particular dynamics that provides a spectral pattern in agreement to data. 
Since existing data beyond 1.9 GeV are plagued with large experimental errors we cannot 
pretend our quantitative predictions to be very precise. Nonetheless we think our findings 
can be of help for the always needed experimental spectral searches in this region as well as 
for the theoretical progress in the dynamical interpretation of data. 
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FIGURES 

FIG. 1. Relative energy nucleon spectrum for the screened potential of Eq. (1) with the param- 
eters of Ref. [7]. The thick solid lines represent our results. The shaded region, whose size stands 
for the experimental uncertainty, represents the experimental data for those states cataloged as 
(* * *) or (* * **) states in the Particle Data Book [13]. Experimental data cataloged as (*) or (**) 
states are shown by short thin solid lines with stars over them and by vertical lines with arrows 
standing for the experimental uncertainties. Finally, we show by a dashed line the Iq ionization 
threshold and by a long thin solid line the total threshold. 

FIG. 2. Same as Fig. 1 for A states. 

FIG. 3. Effective interaction, 3V{rij) — En, for the potentials V{rij) of Eq. (1), dashed line, 
and Eq. (3), solid line, for a two-particle spin 1 state. Ej^ stands for the eigenvalue of the nucleon 
ground state with the corresponding potential. 

FIG. 4. Same as Fig. 1 for the screened potential of Eq. (3) with the parameters of Table I. 

FIG. 5. Same as Fig. 4 for A states. 
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TABLES 



TABLE L Quark model parameters. 



ruu = rud (MeV) 


337 


ro (fm) 


0.495 


Tsat (fm) 


2.12 


K (MeVfm) 


10.0 


Ko- (MeV imj 


1 on n 
izU.O 


a (MeVfm^^) 


976.56 


Mo (MeV) 


- 1726.78 


TABLE IL One-quark ionization thresholds. We 


give the energy above the nucleon ground 


state mass. I is the orbital angular momentum of a 


pair of quarks, s its relative spin and t its 


isospin. 






E (MeV) 


(0,0,0) 


984.10 


(0,1,1) 


1143.33 


(1,0,1) 


1384.15 


(1,1,0) 


1416.48 
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TABLE III. Positive parity N and A ground states for different dominant spatial-spin configu- 
rations up to ~ 3 GcV. The assignment of dominant configurations above 2.4 GcV corresponds to 
an educated guess. Experimental data are from PDG [13]. Stars have been omitted for four-star 
resonances. States denoted by a question mark correspond to predicted resonances that do not 
appear in the PDG (their predicted masses, also indicated by a question mark, appear in Table 
VII). 





S = 1/2 


Model (MeV) 


Exp. (MeV) 


S = 3/2 


(0,0, [3]) 


N(l/2+) 


940 




940 










1232 




1232 


Af3/2+) 


(2, 2, [3]) 




1722 




1680 1720 










1921 




1950 


A(7/2+) 


(4,4, [3]) 


N{9/2+) 


2378 




2220 










2175 




2420 


A(ll/2+) 


(6, 6, [3]) 


iV(13/2+)(**) 






2700 














2950 


A(15/2-)(.*) 


(2,0, [21]) 


A(l/2+) 




1849 




1750 




(2, 2, [21]) 




1938 




1990 




A^(7/2+)(**) 


A(5/2+) 




1871 




1905 




(4,4, [21]) 








? 




iV(ll/2+)(?) 


A(9/2+)(**) 




2193 




2300 




(6, 6, [21]) 








? 




Ar(15/2+)(?) 


A(13/2+)(?) 








? 
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TABLE IV. Negative parity N and A ground states for different dominant spatial-spin config- 
urations up to ~ 3 GeV. The assignment of dominant configurations above 2.4 GeV correspond to 
an educated guess. Experimental data are from PDG [13]. Stars have been omitted for four-star 
resonances. States denoted by a question mark correspond to predicted resonances that do not 
appear in the PDG (their predicted masses, also indicated by a question mark appear in Table 
VII). 



1 J. 1. J J-J ^ W U 1 1 V 1 1 V^V 1 U } 


5 = 1/2 


Model (MeV) 


Exp. (MeV) 


5 = 3/2 


(1,1, [21]) 


A^(3/2-), A^(l/2-) 

V/ /'V/ / 


1410 




1520,1535 








1596 




1675 




iV(5/2-) 


A(3/2-),A(l/2-) 




1517 




1700,1620 




(3,3, [21]) 


N{7/2-) 


2275 




2190 








2153 




2250 




N{9/2-) 


A(7/2-)(*) 




2153 




2200 




(5, 5, [21]) 


iV(ll/2-)(* * *) 






2600 
















A-(13/2-)(?) 


A(ll/2-)(?) 








? 




(3,1, [3]) 






2114 




1930 


A(5/2-)(* * *) 


(5, 3, [3]) 






2153 




2400 


A(9/2-)(**) 


(7,5, [3]) 










2750 


A(13/2-)(**) 



TABLE V. N and A approximate degeneracies up to ~ 3 GeV. Experimental masses (in MeV) 
from PDG [13]. Predicted masses are signaled by a question mark. 



N 


A 


Ar(7/2-)(2190) 


A(7/2-)(2200) 


A^(ll/2-)(2600) 


A(ll/2-)(2650?) 


A^(7/2+)(1990) 


A(7/2+)(1950) 


iV(ll/2+)(2450?) 


A(ll/2+)(2420) 


iV(15/2+)(2900?) 


A(15/2+)(2950) 
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TABLE VI. Ground and first non-radial excitation (denoted by a black dot) correspondence ac- 
cording to our model up to ~ 3 GeV. Experimental masses (in MeV) from PDG [13]. Experimental 
unknown masses are signaled by a question mark. 



N 


A 


Af(5/2+)'(2000) 
Ar(7/2+)'(?) 
7V(9/2+)*(?) 

iV(ll/2+)*(?) 

iV(13/2+)*(?) 


iV(7/2+)(1990) 
Ar(9/2+)(2220) 

iV(ll/2+)(?) 
iV(13/2+)(2700) 

iV(15/2+)(?) 


A(3/2+)*(1920) 
A(5/2+)*(2000) 
A(7/2+)*(2390) 

A(9/2+)*(?) 

A(ll/2+)'(?) 

A(13/2+)'(?) 


A(5/2+)(1905) 
A(7/2+)(1950) 
A(9/2+)(2300) 
A(ll/2+)(2420) 

A(13/2+)(?) 
A(15/2+)(2950) 


A^(3/2-)'(1700) 
iV(5/2-)'(2200) 

Ar(7/2-)'(?) 

Ar(9/2-)'(?) 

A'(ll/2-)*(?) 


iV(5/2-)(1675) 
iV(7/2-)(2190) 
Ar(9/2-)(2250) 
iV(ll/2-)(2600) 

A^(13/2-)C?) 


A(5/2-)'(2350) 
A(7/2-)-(?) 
A(9/2-)-(?) 

A(ll/2-)'(?) 


A(7/2-)(2200) 
A(9/2-)(2400) 
A(ll/2-)(?) 

A(13/2-)(2750) 



TABLE VIL N and A predicted states in the interval [2.2, 3.0] MeV. We denote by a black dot 
the first non-radial excitation. 





N 


A 


J = 


-- 7/2 


A^(7/2+)*(2220) 


A^(7/2- 


)*(2250) 




A(7/2- 


)'(2400) 


J = 


= 9/2 


iV(9/2+)*(2450) 


A^(9/2- 


)*(2600) 


A(9/2+)'(2420) 


A(9/2- 


)'(2650) 


J = 


11/2 


iV(ll/2+)(2450) 








A(ll/2- 


")(2650) 






iV(ll/2+)'(2700) 


Ar(ii/2- 


)'(2650) 


A(ll/2+)'(2850) 


A(ll/2- 


)*(2750) 


J = 


13/2 




A^(13/2- 


-)(2650) 


A(13/2+)(2850) 










iV(13/2+)*(2900) 






A(13/2+)'(2950) 






J = 


15/2 


iV(15/2+)(2900) 













TABLE VIIL Dynamical parity doublets from our analysis up to ~ 3 GeV. Experimental 
masses (in MeV) from PDG [13]. Predicted masses are signaled by a question mark. The black 
dot indicates the first non-radial excitation. 



N 


Ground states 


First non-radial excitation 


Ar(5/2+)(1680) : Ar(5/2-)(1675) 
Ar(9/2+)(2220) : 7V(9/2-)(2250) 
iV(13/2+)(2700) : iV(13/2-)(2650?) 


A^(7/2+)'(2220?) : Ar(7/2-)'(2250?) 
A^(ll/2+)'(2700?) : Ar(ll/2-)*(2650?) 
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